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SUMMARY 

Transient data resulting from approximate step and sinusoidal dis- 
turbances of fuel flow were obtained on a turbojet engine operated at 
one constant simulated altitude and flight condition. The data were 
analyzed to evaluate various procedures for obtaining engine dynamic 
characteristics . 

Three of the more common methods of analysis chosen for study 
gave the following information? (l) Transient analysis yielded the 
basic dynamic characteristics with minimum running time and computation. 
(2) Fourier analysis yielded engine frequency spectra beneficial in 
control synthesis, but was sensitive to instrumentation and approxima- 
tion errors. Computation time was excessive. (3) AC analysis yielded 
the additional high-frequency characteristics of the engine at the 
expense of longer running time and moderate amount of computation. 

With the use of these three methods, accurate engine dynamics 
data could be obtained for use in high gain controls only if the source 
impedances and engine nonlinearities were considered. Also, under 
certain restrictions and approximations, functional relationships 
could be used to derive turbojet engine transfer functions. 


INTRODUCTION 

Dynamical systems are, in general, described and defined by 
differential equations. In many cases, the experimental determination 
of these equations and their utilization in analyses are very awkward 
or difficult. If, as in the case of the turbojet engine, the dynamics 
can be assumed to be essentially linear, alternative more useful repre- 
sentations can be used. 

The Fourier integral theorem is the basic mathematical technique 
which translates the dynamical behavior of linear systems from the time 
domain to the frequency domain and vice versa. In the time domain, the 
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system dynamics axe described by the time response to transient inputs, 
such as a step or a pulse. The frequency domain is characterized by 
sinusoidal response (amplitude and phase) to a sinusoidal input (refer- 
ence 1, ch. 10). 

There is thus a choice as to the experimental techniques to be 
used in investigating turbojet engine dynamics . The use of transient 
inputs gives directly the time domain information which is the most 
direct description of the quality and nature of the system. The use 
of sinusoidal inputs gives directly the frequency response/ and this 
form is often most useful in analyses, such as closed loop control 
system analyses, stability tests, and effects of random inputs (refer- 
ence 2, ch. 4). 

For ideal systems the various experimental techniques, used in 
conjunction with Fourier analysis to convert from one domain to the 
other, would give identical results. If one of the standard inputs Is 
difficult to obtain, that input which is obtained with facility can be 
used. The response to this input can be translated, theoretically, to 
the response to the desired input (reference 3, p. 37, and reference 1, 
ch. 10). Considerations of deviations from linearity, accuracy of 
instrumentation, noise or "hash, " as well as the practical considera- 
tions of l imi tations of running time and of ease of obtaining and pro- 
cessing data, however, separate the various methods as to their 
applicability. 

The object of this report is to provide a detailed comparison of 
three methods of analysis, carried out at the NACA Lewis laboratory, 
in order to make possible a choice of a suitable method of determining 
engine dynamic characteristics. Data resulting from an approximate 
step function input are analyzed by two methods: (l) Functional 

relationships are observed directly from the data which lead to differ- 
ential equations describing the data, and (2) Fourier analysis of the 
data is performed. A third method uses data resulting from sinusoidal 
disturbances to experimentally define the frequency response functions. 
The Fourier transform is used to show the essential equivalence of the 
methods and to compare them. The limitations of each method of analysis 
are discussed with respect to the accuracy of results, ease of obtaining 
and processing data, effect of engine supply sources on data, and 
economical use of engine and test facilities. 


APPARATUS AND PROCEDURE 
Engine Installation 

Engine . - A turbojet engine was mounted on a wing section which 
spanned the test section of the NACA Lewis altitude wind tunnel. The 
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altitude and flight Mach number were simulated in the following manner. 
Air from a climatic control source was applied to the engine 
through a ram pipe. The engine exhausted into the altitude wind tunnel, 
which approached an infinite sink. The wind tunnel static pressure and 
engine inlet total pressure and total temperature were a measure of 
simulated flight speed and altitude. 

The main components of the engine used in this investigation con- 
sisted of an eleven-stage axial-flow compressor, eight through-flow 
combustors, a single-stage turbine, and a fixed -area exhaust nozzle. 

A sketch of this engine is shown in figure 1. 

Fuel system . - The fuel distribution system was of standard design, 
except for the fuel metering valve. This valve was designed to maintain 
a fixed flow rate for each valve position independent of the fuel inlet 
pressure and combustion chamber pressure. The constant flow rate was 
maintained by regulating the pressure drop across the fuel valve with 
a regulatory bypass connected across the fuel valve. 

An independent source of voltage was used to regulate fuel flow 
by feeding the signal into a positional servomotor which was coupled to 
the fuel valve through suitable gear reduction. Because of the 
characteristics of this fuel flow controller, a linear relationship 
existed between the fuel valve position and input signal. 

The DC voltage level determined the engine operating point, and 
a switching of the DC levels produced approximate step changes in 
fuel flow. A variable frequency voltage generator produced a sinusoidal 
voltage of required amplitude and frequency which was superimposed on a 
DC voltage level. In ‘this manner, a sinusoidally varying fuel flow was 
impressed upon the engine about the desired operating level. 


Instrumentation 

All quantities measured were recorded using instrumentation with 
the following characteristics: (l) linear phase shift and flat fre- 

quency response, (2) sensing elements with suitable dynamic response 
and linear electrical outputs, and (3) signal amplification about an 
arbitrary operating point. Those quantities whose characteristics 
were definitely dynamic were recorded on an oscillograph to obtain the 
wave form of the transient, while the initial and final equilibrium 
points were recorded by photographing an instrument panel. 

Instrumentation for transient data . - The signal outputs of the 
various sensing elements were generally too weak to be recorded 
directly and suitable amplifying equipment was therefore necessary. 
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The instrumentation consisted basically of a multichannel recorder, 
associated amplifiers, and sensing element for each channel. The two 
basic recording systems are shown in figure 2 . 

Engine speed was sensed b y a small DC motor coupled directly to 
the compressor shaft. The motor was used as a generator to convert the 
rotational speed of the engine into an electrical signal which was fed 
into a low- pass, "constant k, " filter (reference 4, p. 329) to remove 
commutator ripple and hash. The resultant signal was amplified and 
recorded. A block diagram of this system is shown in figure 2(a). 

Gas pressures were sensed by commercially available transducers 
which incorporated bellows mechanically connected to a bonded strain 
gage bridge having four active gages. The bridge unbalance was fed 
into a strain amplifier and recorded. The sensing elements (pressure 
transducer and connecting tubing) were designed for the desired 
pneumatic behavior (reference 5) . A block diagram of the system is 
shown in figure 2(b). 

Fuel flow was measured with a differential pressure transducer 
which measured the pressure drop across an orifice. For fuel flow 
changes of the order of 1 percent, the pressure drop was considered 
to vary linearly with the fuel flow through the orifice. Fuel valve 
position was converted into an electrical signal by a potentiometer 
mechanically coupled to the fuel valve. 

Instrumentation for steady-state data . - The photopanel recording 
of equilibrium points was accomplished with a conventional instrument 
panel and a camera rigidly mounted to avoid mechanical vibration. 
Pressures were measured by manometers, engine speed by a chronometric 
tachometer, and fuel flow by rotameters. 

Recording technique . - When the engine was in equilibrium at the 
initial conditions, and prior to the proposed transient, steady-state 
information was recorded at the photopanel. Immediately following 
the transient, information at the photopanel was again recorded, thereby 
providing calibration for each trace of the oscillograph. Linearity 
was assumed for the calibration. 

Generally, two methods of calibration were used, depending on the 
type of transient involved. For a transient in which the initial and 
final steady-state values differed, calibration was obtained by using 
values obtained from the steady-state instrumentation. For a transient 
in which the initial and final values were essentially the same, a 
calibration was obtained by recording steady-state values for two 
different engine conditions while maintaining the instrumentation gain 
settings constant. The oscillograph traces were then calibrated using 
the scale gradient thus obtained and the initial steady-state recording. 
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Procedure 

The turbojet engine was operated at a constant altitude of 
15,000 feet and an inlet Mach number of 0.22 for all conditions reported 
herein. Step disturbances in the fuel valve setting were made through- 
out the engine speed range. The magnitude of these disturbances 
resulted in speed changes of 1 percent rated sea-level speed (0.01 N r ) . 

In addition, step disturbances in fuel valve position were introduced 
about an average operating point of 0.94 N r that resulted in speed 

changes of 0.02 N r , 0.04 N r , and 0.08 N r . 

Sinusoidal disturbances in fuel flow were impressed upon the 
engine in the neighborhood of three operating points: 0.88 N r , 0.94 N r , 

and 0.98 At these speeds the amplitude of the fuel flow disturbances 

was adjusted to a value sufficient to cause a change of ±0.01 N r at zero 

frequency for the three operating points investigated. In addition, at 
0.94 N r the amplitude of the sinusoidal fuel flow was Adjusted to pro- 
duce engine speed changes of ±0.02 N r and ±0.04 N r at zero frequency. 

ANALYSIS 

Step function changes are useful test disturbances in determining 
the nature of response. They are easily applied and often approximate 
certain operating conditions. The condition most controlled systems 
are most likely to encounter in operation, however, is one in which 
the input is arbitrarily varied; and the response of this system to 
such a random disturbance may be considered as a linear summation of 
sinusoidal functions over an appropriate frequency range . For the 
latter case, a sinusoidal disturbance of various frequencies experi- 
mentally defines system performance (reference 1, ch. 4). 


Functional Relationships 

When a turbojet engine has a step input disturbance in fuel flow, 
the responses indicate that the engine is a single capacity system. 

Thus, differential equations can be easily derived that fit the experi- 
mental response curves (references 6 and 7). (All symbols used are 
defined in appendix A.) 

These equations are derived in appendix B for a turbojet engine 
with 3 degrees of freedom. The independent variables, or inputs, 
considered are: primary engine fuel flow, reheat fuel flow, and exhaust 

nozzle area. The final equations which describe engine dynamic behavior 
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are in transfer function form, from which it is possible to transform to 
differential equation form, or frequency response form. 

Engine speed . - From the basic assumptions and mathematical develop- 
ment shown in appendix B, the solution of the differential equation 
describing engine speed response to a step disturbance in fuel flow is 


The speed response is exponential in nature between initial and final 
equilibrium points. The speed gain K nw is a measure of sensitivity 

and is the speed change due to a change in fuel flow. The engine time 
constant t is a measure of engine acceleration and is the time for the 
engine to accelerate to 63 percent of the final speed change. A 
theoretical speed response to a step change in fuel flow is shown in 
figure 3(a) . 

An experimental speed response resulting from an approximate step 
change in fuel flow is shown in figure 4. Examination, shows some 
departure from a theoretical first-order response. Although the 
dominant lag in engine speed caused by the. rotor inertia is present, 
additional lags appear at the start of the transient. These lags are 
much shorter than the speed lag and may be attributed to the following 
causes: (l) instrumentation lag, including the lag introduced by the 

speed signal filter; (2) fuel system lag. or the time required for the 
fuel flow to increase to the flow required at the prescribed fuel valve 
setting; and (3) combustion lag, or the time needed for combustion 
equilibrium to be reached. 

The preceding lags appear as a chain of series connected lags 
which can be considered as sequenced systems represented by lumped 
parameters (reference 8, ch. 7). These lags can be considered dead 
time in total effect. Usable results are obtained by replacing the 
chain with a system whose response to a step disturbance is represented 
by dead time plus a first-order response. As shown in figure 4, the 
speed response was extrapolated exponentially to the initial steady- 
state value, and the difference in time between this value and the 
initiation of the transient (time of fuel valve position change) is 
considered dead time. 

Compressor-discharge total pressure . - The response of the 
compressor-discharge total pressure to an approximate step change in 
fuel flow (fig. 4) is seen to be composed of two parts. There is a 
sudden pressure rise due to a fuel flow change at constant speed, 
and an exponential rise due to speed change at constant fuel flow. 

With the use of the basic assumptions and mathematical development in 
appendix B, the transient response of figure 4 can be described by 
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The pressure gain Kp^ w is the total change in pressure due to a fuel 

flow change, and the initial rise ratio d is the ratio of the initial 
pressure rise at constant speed to the total change in pressure. A 
theoretical pressure response to a step change in fuel flow is shown in 
figure 3 (h) . 

As with the speed response, the pressure response was assumed to 
begin at the xime coinciding with the extrapolated initial point of the 
speed response. The Initial rise is assumed to occur at this time. 


Fourier Representation 


Any time excitation or response function may be specified in terms 
of its frequency components according to the met} is of Fourier. Conse- 
quently, it is possible to deal with the dynamics of linear systems 
entirely in the frequency domain (reference 2, ch. 10). 


The useful part of the recordings (fig. 4) where the function is of 
bounded variation, and starts and ends in steady state, can be resolved 
by the Fourier theorem into individual harmonic oscillations. From 
appendix C, the frequency response function is 


r* T 


KG(ico) = 


r>T 


cos at dy(t) - i 


0 


sin a>t dy(t ) 


T 


r\ q 


cos cob dx(t) - i 

10 o 


sin cot dx(t) 


0 


(C5) 


The two types of integral can be evaluated by digital computers for 
each harmonic frequency co in the form 


KG(ico) = 


A(oo) - iB(co) 
C(co) - iD(co) 


(CIO) 


from which expression the amplitude and phase angle can be determined 
in well-known fashion (appendix C). 
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The validity of this method of transposing from the time domain to 
the frequency domain can he checked by the following functions, which 
are typical of engine response data. Given an input function 

, ^ -t/8 

I(t) = 1 + 0.1 e ' 


and an output function 


o(t) = 1 


(t-0.4) 

4 

e 


the theoretical transfer loci, assuming each function to result from a 
step input, can be plotted for both functions (fig. 5), and the transfer 
loci relating the output to the input are shown in figure 6. Also 
plotted are the results of Fourier analysis of both functions for 
24 discrete frequencies . At low frequencies, the approximation is 
practically exact; and at high frequencies the approximation yields an 
error of 20 percent in the phase angle but negligible error in 
magnitude . 


Frequency Response 


If a system is linear or nearly so, the steady-state frequency 
response uniquely describes its characteristics. The application of 
frequency response analysis to physical systems does not have to 
distinguish between lumped or distributed parameters (reference 9, 
pp. 120-123). 


The frequency response functions which best fit the experimental 
data are (reference 1, ch. 8) : for speed 


KG(ico) 


V 

1+Tioo 


and for pressure 


- %c» 

These frequency spectra may "be graphically represented (figs. 7(a) and 
7(h)) hy plotting the logarithm of the output magnitude, and the phase, 
as functions of the logarithm of the frequency. Typical data obtained 
by the sinusoidal variation of the fuel valve of a turbojet engine 
are shown in figure 8. 
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RESULTS 

Several dynamic characteristics of a turbojet engine were decer- 
mined experimentally using two operational inputs. The transient fuel 
flow input was assumed a step input for the responses from which the 
coefficients of the differential equations were determined , and the 
actual fuel flow input was used in Fourier analysis. A sinusoidal input 
was used in directly determining the frequency response. All gains, 
time constants, and frequencies have been made dimensionless using the 
factors given in the symbol list (appendix A) . 

Coefficients of differential equations . - The experimental indicial 
responses obtained by approximate step disturbances in the fuel valve 
setting were analyzed to yield the coefficients of the differential 
equations describing the behavior of engine speed and compressor-discharge 
total pressure (equations (l) and (2)). The speed gain K^,, pressure 

gain K_ w , and initial rise ratio d are shown in figure 9 as functions 
c 

of engine speed. The variation of engine time constant t with engine 
speed is shown in figure 10. The measured time constant, as discussed 
in the analysis section, is an approximation of engine speed lag and 
does not include the higher-order effects caused by additional lags. 

The results of increased magnitude of step disturbances applied 
about an average speed of 0.94 N r are listed in table I. 

Fourier representation . - The particular indicial responses (fig. 4) 
were subjected to Fourier analysis. The transfer loci of pressure, fuel 
flow, and speed for an assumed step disturbance are plotted in figure 11. 
The particular gains are made arbitrary at zero frequency for illustra- 
tive purposes only. 

The frequency response function relating speed to fuel flow is 
obtained using equations (07) and (08) of appendix C, and is graphically 
defined in figure 12. Similarly, the frequency response function relating 
compressor-discharge total pressure to fuel flow is illustrated in 
figure 13. Figure 4 represents indicial responses caused by an increase 
in fuel valve setting resulting in a final speed change of 0.08 N r 

about an average operating point of 0.94 R r . A Fourier analysis of these 
data yields an engine time constant t/t^ of 0.67, and an initial rise 

ratio d of 0.80. The engine time constant was determined by the inter- 
section of the two asymptotes (fig. 12), and the initial rise ratio was 
determined by the asymptote in the high-frequency region (fig. 13) . 

Frequency response . - Although frequency spectra were obtained from 
time functions with the use of the Fourier integral, the frequency spectra 
are directly obtainable using sinusoidal inputs. Typical frequency spectra 
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thus obtained are illustrated in figures 14 and 15. The results of AC 
analysis for three operating points are listed in table II and the effect 
of increased amplitude of disturbance is illustrated in table III. The 
frequency responses of speed to fuel valve setting for the three ampli- 
tudes are illustrated in figure 16. 


DISCUSSION 

Quantitative information characterizing engine dynamic responses 
was obtained by three methods using two common input disturbances. 

The methods are discussed for their relative accuracies, and the com- 
puting and engine times required. 


Transient Analysis 

Although the solutions of differential equations describing 
physical systems give an understanding of its performance, the method 
is difficult to apply to higher-order systems without undue idealiza- 
tion. If the idealized engine can be represented by a single capacity 
system, indicial responses can yield directly the required coefficients 
of the differential equations. If instrumentation, fuel system, com- 
bustion lags, and source impedances are negligible, the dominant lag 
and sensitivities can be determined directly from the traces (fig. 4). 

The assumptions of appendix B are valid only for small inputs. 
Higher-order effects magnified by large inputs must be considered. 

As an example, the speed trace of figure 4 includes what is considered 
to be dead time, lasting from the initiation of the transient to the 
extrapolated point of initial speed response. Another consideration is 
the reactance of the engine with the fuel and air supplies . The 
measured time lags and engine gains include the effect of source 
impedances; and with the use of the equations developed in appendix B, 
idealized time lags and engine gains may be calculated which assume 
zero impedance sources (fully regulated air and fuel supplies). 

A series of indicial responses will quickly bracket the engine 
operating range and yield the sensitivities and lags for each operating 
point tested. Although engine operating time and data reduction time 
are at a minimum, the following limitations must be considered. 

Engine nonlinearity is illustrated in table I in which increased 
magnitudes of disturbance yielded decreased time lags and decreased 
values for the initial rise ratio. Part of the decrease in engine 
time constant can be explained by the engine reactance with the air 
supply source (equation (B7) of appendix B). However, the major part 
of the decrease is attributed to the variation of the engine time 


2321 


2321 


NACA TN 2634 


11 


constant with speed equilibrium, point. For the range above 0.88 N r the 

effective engine time constant is lower at the initiation of the acceler- 
ation, than the effective time constant at the final equilibrium point. 

The initial rise ratio d decreases with increased amplitude dis- 
turbance (table I). Beaction between engine and supply sources would 
cause increased values (equation (B20) of appendix B). For the test 
procedure used, the increasing magnitudes were accomplished about an 
average point 0.94 that is, the larger steps started progressively 

at lower speeds and ended at higher speeds. From figure 9, it is seen 
that the initial rise ratio decreases with decreasing speed. Therefore, 
the initial pressure rise due to fuel flow change occurs at progressively 
lower speeds for large steps, resulting in lower initial pressure rise 
ratios . 

Although it is desired to use disturbance magnitudes that are as 
small as possible to closely approximate the idealized engine, instru- 
mentation sensitivity and error place a lower limit on the usable 
input magnitudes. With instrumentation sufficiently sensitive to record 
engine perturbations and turbulence, the input magnitudes can be reduced 
only to the point at which the responses approach the order of the 
inherent randomness of the engine variables. 


Fourier Analysis 

* The frequency spectra of a turbojet engine are often used as an 

aid in the synthesis of engine controls, and transformation from the 
frequency domain to the time domain is unnecessary (reference 2, ch. 7). 

In the frequency domain, stability and effectiveness can be deduced 
from simple plots in a method which lends itself to the use of analytical 
and experimental data, or to combinations of both. Although step function 
inputs can be used to judge the quality of response, time responses are 
not essential for synthesis. 

In an effort to keep engine operating time at a minimum ^, ari d yet 
produce engine frequency spectra, Fourier analysis was performed on an 
input and the corresponding output. The resultant magnitude and phase 
spectra are shown in figures 11 and 12. Although engine and equipment 
times are at a minimum, computing and analysis times are at a maximum. 

The increasing amplitude ratio in the low-frequency region is 
thought to be a perturbation outside the engine system, and thus not 
> attributable to engine dynamics. As can be seen from figure 11, the 

transfer locus of measured fuel flow to an assumed step input of fuel 
. flow shows the presence of a low-frequency component. Detailed investi- 

gation of the data of figure 4, which were subjected to Fourier analysis, 
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shows that fuel flow reaches a maximum quickly, slowly decreases to a 
minimum, then gradually increases again with time. The frequency 
spectrum relating pressure to fuel flow shows similar behavior. 

The transfer locus of speed to an assumed step disturbance 
(fig. 11 ) shows a departure from a first-order system. Examination 
of figure 4 shows a time delay starting at the initial point of the 
transient. This delay time, or dead time, masks the lag time of the 
speed response * and the engine appears ae a higher-order system. 

However, the frequency spectrum relating speed to actual fuel flow 
(fig. 12) consists of a magnitude spectrum indicating a first-order 
system, and a phase spectrum indicating a higher-order system. Thus, 
the delay time is evidenced as additional phase shift. This indicates 
that the engine is not a minimal phase system (reference 9, pp. 120-123). 
Both the magnitude and phase spectra are necessary to describe engine 
dynamic behavior. For example, unstable combustion during acceleration 
introduces space derivatives causing additional phase shifts. Especially 
for large disturbances of fuel flow, engine behavior departs from that 
of an idealized engine (lumped constant, minimal phase); and engine 
behavior is described by a distributed parameter system which is not of 
minimal phase. 

Comparison of figures 12 and 14 with figure 7 (a) indicates that the 
speed response to actual fuel flow is exponential in nature and that the 
engine is primarily a first-order system. The comparison of figure 13 
and 15 with figure 7(b) shows that the pressure response to actual fuel 
flow is of a lag nature (reference 1, p. 265). The high-frequency 
scatter in figures 12 and 13 is due to the nature of the approximation 
used (appendix B) in the Fourier analysis. An increase in the number 
of ordinates for approximating the transient recordings would have 
reduced the high-frequency scatter. As shown by figures 5 and 6 this 
particular procedure for Fourier analysis is valid for functions that 
are continuous and differentiable for t a 0 (appendix C); however, it 
is very sensitive to perturbation and noise contained in the actual 
responses, and to the number of ordinates used. The number of ordinates 
per second of transient time must be kept as high as practical. 


Frequency Response 

Because frequency spectra are useful in visualizing the dynamic 
characteristics of a physical system, sinusoidal disturbances were 
applied to the turbojet engine. The frequency of oscillation was 
increased until readability of data waB limited by instrumentation. 
Figures 14 and 15 are plots of the frequency response of speed, pressure, 
and fuel flow to fuel valve position. 
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The magnitude K |G(ico)| is composed of a frequency invariant part 
K and a frequency variant part iG-(ico) { . For this investigation, K 
is the engine gain, relating the fuel valve position with the speed, 
pressure, and fuel flow; G-(ito) represents that part which varies with 
frequency and it is a function of the energy storage of the engine. 
Figures 14 and 15 have been normalized by dividing the magnitude by K 
and therefore apply only to the frequency variant portions of the magni- 
tude. 


In figure 14 the first comer frequency of the speed response 
curve characterizes the engine time constant T. The 45 ° slope defines 
a first-order system, and the second corner frequency is due to the 
lag in the fuel system. Information contained in the higher-frequency 
range demonstrates the nature of the fuel system. If this response is 
approximated by the drawing of a 45° slope through most of the data 
points, an equivalent time constant results. It is seen that the 

fuel system source impedance is reflected in the speed response. 

Figure 15 shows a similar effect. 

The effect of amplitude of disturbance is illustrated by figure 16. 
Increased amplitude yields decreased values of engine time constants. 
Although the magnitude spectrum indicates a lumped constant system, the 
phase spectrum indicates additional phase shift which can result from 
dead time. Thus, the engine is not a minimal phase system. Because 
the amplitude ratio increases with amplitude of disturbance at a 
given frequency, it is possible that there exists a dead band which 
would account for this effect for, in the extreme case, a disturbance 
amplitude entirely within the dead band would result in zero output 
amplitude. 

Engine frequency response appears to yield the same information as 
transient analysis and Fourier analysis, especially in the low-frequency 
domain. The method of frequency response yields additional high- 
frequency data, such as fuel system lag, within the limitation of the 
type of instrumentation used. The nature of nonlinearity such as dead 
time is also revealed. Engine and equipment operating times, however, 
are at a maximum. 

The effect of excessive lags in a control loop for an engine may 
result in instability if the loop gain is too high. Control systems 
become more oscillatory as the loop gain Is increased, and the oscilla- 
tion can be sustained indefinitely or increased with time. Although 
stability criteria may indicate for idealized second-order or first- 
order systems that gain can be increased indefinitely without instability, 
in physical systems this is impossible. Instability will always result 
because of factors neglected in the idealization. If Nyquist's 
criterion did indicate the presence of roots with positive real parts 
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(reference 1, p. 175) , the oscillations resulting from, a small disturbance 
would not, in a physical system, increase indefinitely in magnitude. 

Under these circumstances, the oscillations will increase in magnitude 
until the apparatus destroys itself, or until nonlinearities limit the 
oscillations to some finite magnitude. 

\ 

A recording of controlled engine acceleration using a high gain 
controller in a speed - fuel-flow loop is shown in figure 17. The block 
diagram of the control loop is shown in figure 18. Fuel system lag, 
coupled with combustion delay time and instrumentation lag, was suffic- 
ient to reduce the gain margin to approximately 26 decibels (fig. 14). 

In figure 17 speed is controlled by fuel flow, with an upper limit on 
the fuel valve opening to prevent compressor stall and surge. As top 
speed is approached, the control gain is suddenly increased in an 
effort to increase control effectiveness. The resultant loop gain, in 
the vicinity of 26 decibels, was sufficient to cause oscillatory fuel 
flow at a frequency in the spectral region of the fuel system time 
constant indicated in figure 14. 


Comparison of Methods 

The time constant which resulted from frequency response analysis 
was consistently higher than the time constant which resulted from 
transient analysis. Because of the approximation used in transient 
analysis, the time constant did not include higher-order effects; 
whereas, the apparent time constant which resulted from frequency 
response analysis was affected by the nearness of the equivalent fuel / 

system time constant (fig. 14) . Similarly the values for initial rise 
ratio were decreased values because of the nearness of the equivalent 
fuel system time constant (figs. 14 and 15). A comparison is made in 
table IV. 

Increased amplitude disturbances yielded results which demonstrated 
engine nonlinearity (table III). For increased magnitudes of distur- 
bances, the frequency response method again yielded values of engine 
time constant that were higher than those obtained from transient 
analysis (table I). Thus, both methods of analysis showed similar 
nonlinear effects for evaluating the engine time constant. However, 
both methods of analysis showed dissimilar nonlinear effects for 
evaluating the initial pressure rise ratio. A comparison is made in 
table V. 

Fourier analysis resulted in values of engine dynamic character- 
istics comparable with those obtained by AC analysis but higher than 
those obtained from transient analysis. As explained previously, 
transient analysis used an approximate method which did not include t 

higher-order effects. A comparison is made in table VI. 
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A comparison of the three methods for determining the speed response 
to fuel flow disturbances in the frequency domain is shown in figure 19. 


SUMMARY Of RESULTS 

The three methods of. analysis applied to turbojet engines produced 
engine dynamics data usable to varying degrees in control synthesis. 

The choice of method is dependent on the engine and calculating times 
available, degree of accuracy desired, type of instrumentation used, 
and the application of the results. 

The following results were obtained from the evaluation of the 
analyses : 

1. Transient analysis yielded the basic engine time constant and 
various engine gains with minimum engine, computation, and equipment 
operating time. Under the limitation of the instrumentation and 
recording system, however, transient analysis did not yield the rapid 
response or high-frequency characteristics of the engine and its 
associated systems . 

2. Fourier analysis of engine transient data produced engine fre- 
quency spectra beneficial in control synthesis. However, this method 
was sensitive to instrumentation and approximation errors. Calculating 
time was excessive. It would be especially excessive if the approx- 
imation accuracy were increased in an effort to obtain high-frequency 
data. 

3. AC analysis yielded high-frequency information, such as fuel 
system lags, at the expense of engine and equipment operating time 
and of calculating time. This type of information was necessary in 
high-speed and high gain control systems. 

4. Only if the source impedance of the test facilities and the 
nonlinearities of the turbojet engine were considered could accurate 
engine dynamics data be obtained for use in high gain controls. 

5. Under certain restrictions and approximations, functional 
relationships could be used to derive turbojet engine transfer func- 
tions which fit the experimental data and which are easily manipulated 
in control loop synthesis. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 26, 1951 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

A exhaust nozzle area, sq. ft 

A,B,C,D constants 

d initial pressure rise ratio due to step change in fuel flow 


e 

f 


g 


G 

I 

J 

E 



at constant A, Wg, and P; — 


N 


P c w 


base of natural logarithms 
initial pressure rise ratio due to step change in reheat fuel 

r dP, 


flow at constant A, W, and P; 


(—) 


initial pressure rise ratio due to a step change in area at 
constant W, Wg, and P; 



N 


Sc* 


initial pressure rise ratio due to a step change in inlet total 

/a p„\ 


l 


ap 


pressure at constant W, and A; 

functional relationship 
input function 

polar moment of inertia, (lh-ft) (sec^) 
frequency invariant gain 


N 


*p, 


i 
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■^na 

speed 

dN 

dA' 

^np 

speed 

W 2 ' 

K nw 

speed 

dN 

dW' 

^2 

speed 

and 

K p c a 

pres si 
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E 


■PfiP 




E 


Pc w 2 


speed gain due to area change at constant W, Wg, and PJ 
, rpm/sq. ft 


pi ij lh /hr 

a in due to reheat 


dP c dN f SPc> 


dA dA 


lbM ft / B5 ft 


pressure gain due to inlet total pressure change at constant W, 

^c _ dN (&c\ fdPc\ 

W 2 > and A; gp dP J ^ J r 


N 


pressure gain due to fuel flow changes at constant Wg, A, and P; 

w 2 - §(irl + (I 2 ),,’ lb/sl f /“ /hr 


pressure gain due to reheat fuel flow changes at constant W, 

dP c _ dN 
A, and P; dW d u 

^ Wo ^ 'N 



+ (*o\ 

Van I 

Uw A 


, lb/sq. ft/lb/hr 
im/lb /hr 

(K n _J 0 . 94 N r Pressure gain evaluated at 0.94 N r , lb/sq. fy^Lb/hr 

, lb/sq f^i 


(K nw )o 9 speed gain evaluated at 0.90 N r , rp: ^ t 


Pc w 
k 

N 


proportionality factor, 
Laplace transformation 


'rpm 


engine speedy rpm 


I 
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*r 

0 

P 

P 

Q 

3 

T 

t 

W 

Wc 


rated sea-level engine speedy rpm 
output function 

engine inlet total pressure, lb/sq ft 
compressor-discharge total pressure, lb/sq ft 

number of cycles in time T, cycles 

engine torque, ft-lb 

Laplacian operator 

total time, sec 

instantaneous time, sec 

engine primary fuel flow, lb /hr 

engine reheat fuel flow, lb/hr 


X 

input function 


Y 

variable 


y • 

output function 


A 

incremental change 


CT 

correction factor, ^1 + k 

dN\ 

% 

engine time constant, sec 



engine time constant obtained by AC analysis for engine 
speed of 0.88 N r , and amplitude of +0.01 N r , sec 

equivalent fuel system time constant, sec 


<P 

co 


phase angle, radians 
circular frequency, radians/sec 
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APPENDIX B 

DERIVATION OF TRANSFER FUNCTIONS DESCRIBING 
ENGINE DYNAMIC BEHAVIOR 

Examination of transient data suggests that engine Behavior can he 
described by two general types of transfer functions of the first order 
(reference 2, ch. 5), the dominant characteristic of which is a speed 
lag. The simultaneous relationships of pressure, temperature, fuel 
flow, and air flow lead to quasi-static assumptions, at least for small 
deviations from the equilibrium point (references 6 and 7). If certain 
functional relationships are assumed, simple transfer functions result 
that agree with experimental data. These transfer functions are 
restricted by the following assumptions: 

(1) Engine is considered a lumped constant system. 

(2) Linearization about an arbitrary operating point applies. 

(o) Conditions of flight speed and altitude are constant. 

(4) Only major storage element is rotor inertia, J. 

(5) All thermodynamic and flow processes are quasi-static. 

(6) All supply sources have zero impedance. 

Engine Speed 

The experimental speed response of figure 4 to an approximate step 
disturbance in fuel flow indicates an exponential relationship. If the 
assumption is made that 


and that linearity applies 
equilibrium point, then 

AQ 


For a turbojet engine 
increase in fuel flow produces an unbalance torque which accelerates 
the engine rotor to a new equilibrium speed. Using the relationship 


Q = G(N,W) (Bl) 

for small deviations from a steady state or 



initially in steady state, an incremental 
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AQ = JAN 

and equating (B3) with (B2) yield 


(B3) 


JAN 


= (— ) AW + (~) AN = J ~ 

V3w 4 at 


AN 


or 


( 


jfc\ 


AN = - 


W 


■AW 


(B4) 


1 - 


at 

Un/ w 


Replacing rr hy the Laplacian operator, for the Initial condition of 
at 

steady state, equation (B4) Becomes, in transfer function form (refer- 
ence 8, ch. 5), 

^nw 


AN(s) = 


1 + Ts 


AW(s) 


(B5) 


where the speed gain K nw is ^ and the engine time constant T 

-J 


is 


( 


oN/ 


W 


By definition, the transfer function is the ratio of the Laplace 
transform of the output to the Laplace transform of the input; or from 
equation (B5) 


AN / ■> ^nw 
~ aw S " 1 +Tb 


(B6) 


If the output results from a step input disturbance, then the responses 
are called indicial responses. The excitation function, or step distur^ 
Bance, multiplied By the transfer function (equation (B6)) results in a 
response transform 


AN(s) 


AW ^nw 

S 1 + TS 


whose inverse transformation is 
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(B7) 

A theoretical indicial speed response is shown in figure 3, which compares 
with the experimental indicial response of figure 4. The indicial response 
function completely defines, as well as the transfer function, the relations 
of a physical system. If the indicial response is known, it is possible 
to compute with its help the effect of any arbitrary disturbance which 
affects the system. 

This speed - fuel-flow relationship applies for an idealized engine 
with only one degree of freedom (constant exhaust nozzle area, zero reheat 
fuel flow, and constant inlet total pressure) and subject to the listed 
restrictions . 

Obviously, for changes in only one independent variable, with the 
others held constant, the following transfer functions apply: for area 

changes 


AN(t) = 


(i - .■*/*) 


M. ( a ) = £ria 

AA 1 + Tg 

for reheat fuel flow changes 

AW / v ^ TlV 2 

AW 2 (BJ 1+^8 

and for inlet total pressure changes 


(B8) 


(B9) 


M ( S ) _ 

AP v ' " 1 + Ts 


(BIO) 


It follows from the linearization assumption that the combined speed 
response resulting from simultaneous variations of input variables is 
the superposition of. responses from variations in each independent 
variable separately. Thus the addition of equations (B6), (B8), 

(B9), and (BIO) is the combined speed response, or 


AN(s) 


3™ . AW(s) + ■ Kn I 2 - AW 2 (a) 

1 + TS 1 + TS 


+ — 3 a ~- AA(s) 

1 + TS 


+ AP (s) 

1 + Ts 


(Bll) 


where the E's are defined in appendix A. 

It also follows from the linearization assumption that the engine 
time constant x is the same regardless of which independent variable 
produces acceleration. 
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The idealized transfer function (equation (Bll)) assumes zero 
impedance sources; that is, there must he no reaction between the engine 
and the independent variables. However, in turbojet engine test instal- 
lations, 100 percent regulation of fuel and air supplies is difficult. 
Fuel flow does not respond instantly to fuel valve setting because of 
the necessary plumbing; and the variation in combustion chamber pressure 
during engine acceleration varies the pressure drop across the fuel 
nozzles and hence varies the fuel flow. 

Similarly, air is supplied through large ducting and mass flow is 
adjusted by throttling. Regulation during acceleration is difficult 
and engine reactance with air supply is noticeable as indicated by 
figure 4. If the change in inlet total pressure is assumed proportional 
to speed. 


AP = - kAN (B12) 

The proportionality constant (-k) is negative if the pressure P 
decreases for increases in engine speed. For constant exhaust nozzle 
area and zero reheat fuel flow, equation (Bll) becomes 

“W - t'Vto ™ <b) - rvfe kA1,<B) 

or 


AN 

AW 


(a) = 


E; 


■nw 


1 + 


^p 


1 + 


^TLt 



(B13) 


Thus, the experimental indicial response (fig. 4) includes the variation, 
during acceleration, of the throttling losses on the air supply. Com- 
parison of equation (B13) with equation (B6) indicates that the time 

constant - and speed gain are reduced values because a = 1 + hE™ 

a C7 

is always greater than one. 


The transfer function yields the frequency response function with 
the substitution s = io> (reference 8, ch. 3). If the derivatives are 

substituted for the E's and for s, the transfer function 

dt 

becomes a linear differential equation. 
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Compressor-Discharge Total Pressure 


The experimental pressure response of figure 4 indicates that there 
is an initial pressure rise at constant speed due to a fuel flow change, 
and a subsequent pressure rise at constant fuel flow due to the resultant 
speed change. If, 


P c = G(N,W) 


(B14) 


and linearity is assumed for small excursions from the steady state 
point, then 


AP, 


f& c \ [i%>\ 

= — -) AW + — £ AN 

VcJW /jj \ dN / w 


From equation (B4) 


APf 


V 5 W J 


AW + 


a_N ft P c\ dN f dF c\ + c\ (ftp) d 

^W\Sn/ w dW \bs VcW VdW / dt 

AW = :: 


AW 


1 +T dt 


1 + T it 


« K, 


P_» 


1 + dT st 
1 +T it 


AW 


(B15) 


wtL©r© E- v -I- ( = 

awv^iM ) u Vow ; N 


dN c | 

3 C W ■ <Mbn i w 

r 


dPc 

dW 


= pressure gain 


and d = 


N 


E 


= initial pressure rise ratio. In transfer function 


p n w 


form, equation (B15) becomes 

AP 


( \ - 1 + dTs 
AW l ' 3 ' 1 + ts 


(B16) 


The response transform to a unit step excitation function is 


AP c (s) = 


1 + dxs 1 

1 +TS 3 




whose inverse transformation is 
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AP c (t) = Kp cW AW [jL - (1 - d)©"" 1 / 17 ] (B18) 

A theoretical indicial pressure response which compares with the experi- 
mental indicial response of figure 4 is shown in figure 3. 

Again, as with the combined speed response, it follows from the 
linearization assumption that the pressure response resulting from 
simultaneous variations of input variables is 

“e<«> - %„« {Hf?} “ <3) + {Hvr} AW ^ (s) + 

K p c a {Hlr} + *s c p {HI?} 4P < 3 > (B19) 

The values of the K‘s and d, f, g, and h are defined in appendix A. 

If, again, the change in inlet total pressure is assumed to be pro- 
portional to the speed change, 



or, with the use of equation (B13), 



o 


(B20) 

where the experimental pressure response now includes the variation of 
the throttling losses on the air supply, or 
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Wp ' k (& C \,tL 


N,P 


In comparison with equation (B16), the initial rise 



remains 


the same "because it is the pressure rise due to fuel flow changes and 
it occurs initially at constant speed. The additional pressure rise. 


of equation (B16) . This additional pressure rise, at constant fuel 
flow, due to speed change is a reduced value. Consequently the initial 


is an increased value. 

The observed responses of any pressure, temperature, or thrust 
are similar to the response of compressor-discharge total pressure. 
Thus, the development of the respective transfer functions is similar 
to the foregoing development. 



pressure rise ratio 



d 
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APPENDIX C 

FOURIER ANALYSIS OF TRANSIENT DATA 


It is possible to transform, from the time domain to the frequency 
domain if an arbitrary input and corresponding output are known. If 
the input and output are related by a transfer function, then 


KG ( S ) = = & [output] 

-S^[x(t)] y [input] 

For a system initially at rest, 


■?[ y(t)] = sy[y(t)] 


so that in this case 


KG(s) 


sif Cy(t)] 

s y[x(t)] 


y(t)e- st 


x(t )e -3 ^ 


dt 


dt 


(Cl) 


(C2 ) 


where the numerator and denominator are unilateral Laplace transforms 
of y(t) and x(t). 


Because e -3 ^ is continuous and y(t) is of bounded variation in 
(0,»), the Stieltjes integral of e" st with respect to y(t) from 0 
to •» exists (reference 10, p. 7). 


If both the numerator and denominator of equation (C2) converge, 
they define a function of s called the Laplace-Stielt jes transforms 
of y(t) and x(t), respectively (for example, reference 10, ch. 2). 


G-(s) 


m 



d[y(t)] 


( 03 ) 


However, if 


G(s) = 


e" 3t y(t)dt 


then G(s) is the Laplace transform of y(t). 
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If y (t ) and x(t ) exist for t>0 and approach, zero sufficiently 
fast such that y(t) and x(t) are approximately constant for t>T, 
then equation (Cl) becomes 


KG(s) = 


J ^T 

0 


e -£rb y(t)dt 


'T 


(C4) 


e“ st x(t)dt 


By the transformation s = ico, equation (C4) can be written as a 
Fourier-Stielt jes transform (reference 10, p. 252); and the frequency 
response function can now he expressed as (reference 11 ) 


KG(ico) = 


r--*' 


y(t)dt 


rv-* 


x(t)dt 


P T . 1 

PT 

1 y(t)cos cot dt - i 1 

y(t)sin cot dt 

Jo 

10 

pT 


1 x(t)cos cot dt - i 

x(t)sin cot dt 

Jo J 

0 


Thus there are two types of integrals to be evaluated 

Re[y(tj] - i Im[y(t)] 


KG(ico) = 


ReJjx(t)] - i Im[x(t}f] 


(C5) 


(C6) 


Consider 


Let 



rvji 

Be[y] » 

y(t) cos cot dt = 

j 

0 


ii 

£ 

i 


r»T 


0 


cos cot dy 


where y^ is the final value, then 


dY = 0 - dy 



At t = 0, 
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at t = T, 


y = y f 

Y = y f - y f = o 


Thus 


H>] = - f 


0 


r\ T 


cos ®t dY = 

T w 


Integration by parts gives 

Be[y] = Y cos cot 


cos cot dY 


0 


+ co I Y sin cot dt 


= Y 0 “ “ 


'T 


Y sin cot dt 


(C7) 


Similarly, 


nia[y] = Y sin cot 



up 


(C8) 


Let T = total transient time and p = number of cycles in time T. 
The domain 0 $.t $T is divided into 120 equal intervals; and the 
integration is done by using a seven-point numerical integration 
scheme (reference 12, pt. II, p. 233), in conjunction with a digital 
computor, 



ife [ 4 


41y 0 + 216y 1 + 27y 2 + 272y 3 + 27 y 4 + 216y 5 + 41y ( 


0 


(C9) 


then 


radians x p cycles ^ 12Q ordinates _ 2itp radians 
cycle 120 ordinates T sec T sec 
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and 


or 


and 


where 


Re[y] = Y, 


l 0 


- 2-iffi A t sin S? t 

T 120 z 60 

t=0 


120 


He 


120 

t=o 

jtp . 
60 t 

120 



^0 = lo ^ AtCOS ^0 


. 41 41 , > 

*0 = 335 Y ° ■ uo (y f • y o> 


. 216 216 , v 

A i ■ ii5 Yl * H5 (yf ' yi) 


. 27 27 , . 

A 2 = iio Ys * ns (yf ' y 2 


Thus the frequency response function can he readily evaluated for each 
frequency o>. 

- cfei : Si“i < oio > 

The amplitude and phase are determined in -well-known fashion as 


K 


|g(1co)| =J 


A 2 -h B 2 
C 2 + D 2 


(Cll) 


<p = tan' 


■1 AD - BC 
AC + BD 


(C12) 
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TABLE I - VARIATION OF ENGINE DYNAMIC 
CHARACTERISTICS WITH MAGNITUDE 
OF DISTURBANCE 

[Average operating point =0.94 N r ] 



Magnitude of 
disturbance 
(percent N ) 

Engine time 
constant 

t/t i 

Initial pressure 
rise ratio, d 

1 

0.78 

0.80 

2 

.74 

.80 

4 

.68 

.75 

8 

.61 

.71 


TABLE II - VARIATION OF ENGINE DYNAMIC 
CHARACTERISTICS AS DETERMINED 
BY AC ANALYSIS 

[Amplitude of disturbance at zero frequency = ±0.01 N r ] 


Operating point 
(percent N r ) 

Engine time 
constant 
t 

Initial pressure 
rise ratio, d 

98 

0.67 

0.81 

94 

.87 

.70 

88 

1.00 

.54 
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TABLE III - VARIATION OF ENGINE DYNAMIC 
CHARACTERISTICS WITH MAGNITUDE 
OF DISTURBANCE AS DETERMINED 
BY FREQUENCY RESPONSE 


[Operating point =0.94 N r ] 



Amplitude of 
disturbance 
(percent N^) 

Engine time 
constant 

t/t 1 

Initial pressure 
rise ratio, d 

±i 

0.87 

0.70 

±2 

.86 

.72 

±3.5 

.72 

.78 

±6 

.67 

.80 


TABLE IV - COMPARISON OF ENGINE DYNAMIC 
CHARACTERISTICS OBTAINED BY TWO 
METHODS OF ANALYSIS FOR 


0.01 N DISTURBANCES 
r 


Transient analysis 


AC analysis 

Operating point 

Engine time 

Initial 

Engine time 

Initial 

(percent N r ) 

constant 

rise 

constant 

rise 

i — i 

ratio 

d 

i — i 

ratio 

d 

98 

0.61 

0.92 

0.67 

0.81 

94 

.78 

.80 

.87 

.70 

88 

.98 

.67 

1.00 

.54 
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TABLE V - COMPARISON OF ENGINE DYNAMIC CHARACTERISTICS 


OBTAINED BY TWO METHODS OF ANALYSIS FOE VARYING 
MAGNITUDES OF DISTURBANCES 


[Average operating point = 0.94 N r ] 



Transient analysis’ 

AC analysis 


Magnitude of 

Engine time 

Initial 

Amplitude of 

Engine time 

Initial 

disturbance 

constant 

rise 

disturbance 

constant 

rise 

(percent N r ) 

t/t 1 

ratio 

d 

(percent N r ) 

i — 1 

e 

ratio 

d 

1 

0.78 

0.80 

±1 

0.87 

0.70 

2 

■ .74 

.80 

±2 

.86 

.72 

4 

.68 

.75 

±3.5 

.72 

.78 

8 

.61 

.71 

±6 

.67 

.80 


TABLE VI - COMPARISON OF ENGINE DYNAMIC CHARACTERISTICS 
OBTAINED BY THREE METHODS OF ANALYSIS 
[Average operating point =0.94 N r ; 

approximate magnitude of disturbance = 0.08 N r J 


Method of analysis 

Engine time 
constant 
t/t-l 

Initial 

rise 

ratio 

d 

Transient 

0.61 

0.71 

Fourier 

.67 

.80 

Frequency response 

.67 

.80 


0 1 
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su pply/ 



Combustion 
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Turbine 


Fixed-area 
exhaust nozzle 



Exhaust into 
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wind tunnel 


1 w 

Fuel 

valve 

discharge 

total pressure 

1 

P c 

Fuel 

supply 


w 


Proportional 

servomotor 


Input signal 
(volts) 


Figure 1. - Turbojet engine installation illustrating supply sources and control, 

and location of sensing elements. 
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N(rpm) 


P(in. Hg) 


Figure 2. 


Speed sensor Filter 


Amplifier and recorder 



Maximum pen 
sensitivity 



(a) Engine speed. 


Pressure sensor 


Pneumatic 

system 

1 

Transducer 
din. 
25 inT 

+ 1,2 S+] 
2500 50 b 

IV in. Hg 


Strain amplifier and recorder 


n 0.4 mm pen displacement 
Hin./in. strain 


DC to 100 cps 
(essentially flat) 




0.1 


in. Hg 
mm 


(b) Pressure. 

- Dynamic characteristics of two basic sensing systems for recording transient engine behavior. 
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(a) Engine speed. 



(b) Compressor-discharge total pressure. 

Figure 3. - Theoretical indicial responses of pressure and 
engine speed to fuel flow. 
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mation by digital com- 
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Figure 7. - Frequency spectra of a first-order dynamic system. 
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CT> 



(b) G(s) = L . t 

1 + s 


Figure 7. - Concluded. Frequency spectra of a first-order dynamic system. 
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Figure 9. - Variation of several engine gains with speed for a turbojet engine operating 
at simulated altitude of 15,000 feet and simulated flight Mach number of 0.22. 
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Dimensionless time constant 





Dimensionless speed, ■ — 

N r 


Figure 10. - Variation of engine time constant with speed for a turbojet 
engine operating at simulated altitude of 15,000 feet and simulated 
flight Mach number of 0.22. 
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(a) Amplitude ratio. 


Figure 12. - Frequency spectrum of engine speed response to fuel flow, obtained by Fourier 
analysis of transient responses of figure 4, of turbojet engine for simulated altitude 
of 15,000 feet and simulated flight Mach number of 0.22. 
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(b) Phase shift. 


Figure 12. - Concluded. Frequency spectrum of engine speed response to fuel flow, obtained by Fourier 
analysis of transient responses of figure 4, of turbojet engine for simulated altitude of 15,000 feet 
and simulated flight Mach number of 0.22. 
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(a) Amplitude ratio. 



(b) Phase shift. 


Figure 13. - Frequency spectrum of compressor-discharge total pressure response to fuel flow, obtained 
by Fourier analysis of transient responses of figure 4, of turbojet engine at simulated altitude 
of 15,000 feet and simulated flight Ma f ch number of 0.22. 
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Figure 14. - Frequency spectra of fuel flow and engine speed responses with fuel valve position 
varying sinusoidally at constant amplitude for turbojet engine operated at simulated altitude 
of 15,000 feet and simulated flight Mach number of 0.22. 
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Figure 15. - Frequency spectrum of compressor- discharge total pressure response with fuel valve position 

varying sinusoidally at constant amplitude for turbojet engine operated at simulated altitude of 15,000 feet 
and simulated flight Mach number of 0.22. 
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Figure 16. - Engine speed frequency response to fuel valve position varying sinusoidally at three amplitudes 
for turbojet engine operated at simulated altitude of 15,000 feet and simulated flight Mach number of 0.22. 
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(b) Phase shift. 


Figure 16. - Concluded. Engine speed frequency response to fuel valve position varying sinusoidally at three 
amplitudes for turbojet engine operated at simulated altitude of 15,000 feet and simulated flight Mach 
number of 0.22. 
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Figure 18- - Block diagram of engine-speed - fuel-flow control loop including limits. 
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Figure 19. - Comparison in frequency domain of three methods for determining engine speed response to fuel 
valve position disturbance for a turbojet engine operating at 0.94 N r , simulated altitude of 15,000 feet 
and simulated flight Mach number of 0.22. y 
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Figure 19. - Concluded. Comparison in frequency domain of three methods for determining engine speed response to 
fuel valve position disturbance for a turbojet engine operating at 0.94 N r , simulated altitude of 15,000 feet, 
and simulated flight Mach number of 0.22. 
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